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The internal vortex structure of a trapped spin-1 Bose- 
Einstein condensate is investigated. It is shown that it has a 
variety of configurations depending on, in particular, the ratio 
of the relevant scattering lengths and the total magnetization. 
PACS number: 03.75Fi 

Recently the MIT group has succeeded in obtaining 
Bose Einstein Condensation (BEC) of 23 Na atoms in 
an optical trap, pig] A novel aspect of this system is 
that 23 Na atoms possess a hyperfine spin, with / = 1 
in the lower multiplet. All three possible projections 
of the hyperfine spin can be optically trapped simulta- 
neously. Thus generally the condensate have to be de- 
scribed by a spin-1 order parameter consisting of the 
spinor (\&„, ty , ^fa) where \fr u ,o,d ar e the macroscopic 
wavefunctions with the hyperfine projection s = 1,0, —1 
respectively. The ground state of the condensate is deter- 
mined by the spin-dependence of the interaction between 
the bosons. |3|,|j In the dilute limit they can be charac- 
terized by the (s-wave) scattering lengths ciq and a-i in 
the total (hyperfine-) spin and 2 channels respectively. 
The polar state, where the order parameter is propor- 
tional to the spinor (0, 1, 0), is favored if a 2 > ao, but the 
axial state, where the spinor is (1, 0, 0), if the inequality 
is otherwise. 23 Na belongs to the former case. However, 
87 Rb is predicted to belong to the latter case. || 

A particular interesting feature of these systems is 
that, due to the weak spin dependence of the interatomic 
interaction, ao ~ a 2 - For example for 23 Na ao « 46a^ 
and a 2 rs 52as, where as is the Bohr radius. Thus the 
difference a 2 — a$ rs 6as is only a small fraction of ao 
or a,2 . As a result, the condensate only needs to pay 
very little extra energy to get into the "wrong" state. If 
there is a competing energy, such as that due to the pres- 
ence of a gradient, it may be energetically favorable for 
the condensate to deviate locally from the polar state. 
An analogous remark is also applicable to 87 Rb, where 
ao ~ llOas and a 2 « 107as- [|5| 

In this paper I shall illustrate this by considering vor- 
tices of this spin-1 Bose condensate. Vortices with a 
scalar order parameter in BEC have been discussed in 
many papers (e.g [pP] ). In a cylindrically symmetric 
trap, there is no stable vortex state for angular momen- 
tum L < Nh where N is the number of particles. A 
singly quantized vortex has L — Nh, with the node of 
the order parameter located at the center of the trap. 
Here I shall show that vortices in a spinor condensate 
are even more interesting in that they exhibit a very rich 



internal structure. In general locally the order parameter 
is in neither the polar nor the axial state. They may have 
broken cylindrical symmetry with nodes of the order pa- 
rameter of individual species appearing at positions other 
than the trap center. The minimum angular momentum 
required for the formation of a vortex is also less than 
N. Moreover transitions between different internal vor- 
tex structures are possible. Vortices were also discussed 
in Ref || and 0, but they did not consider structures 
which arise from deviations of the order parameter from 
the original polar or axial phases. 

Since the order parameter has more than one compo- 
nent, it is convenient to distinguish between vortices of 
individual order parameter component $s and the com- 
posite structure. I shall refer the latter as the composite 
vortex (CV). 

The order parameter, in particular that of the CV, 
is found by minimization of the energy (restricting our- 
selves to T = 0) under appropriate constraints. The en- 
ergy density E consists of the kinetic and potential con- 
tributions J2 S ' 2M S ' — r- Vl^sl 2 , where M a is the atomic 
mass, V is the trap potential and the sum is over all spin 
components, and the interaction part which can be writ- 
ten as £i„ t = i(co + c 2 )n 2 - \c 2 \2^ u ^ d - * | 2 where 
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n = E s |*« 
and C2 = — 



is the local density. Here cq = 



—_ 90+292 
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.. where <?o,2 are in turn related to the 
scattering lengths in the total spin and 2 channels via 
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The total particle number TV, angular 
momentum L and magnetization M should be considered 
as conserved if no exchange of the corresponding quantity 
is allowed between the atoms inside the trap and their en- 
vironment (within the relevant experimental time scale) . 
In this case we have to minimize the total energy for given 
N, L and M. As usual it is convenient to introduce and 
minimize the free energy G = E — jiN — HL — HM where 
H, fl and H are Lagrange multipliers, /i, Q, H correspond 
to the chemical potential, angular velocity and magnetic 
field. f2 is given by the angular velocity of the rotating 
trap if angular momentum can be exchanged between the 
trapped atoms and their environment. 

It is useful to note that the energy is invariant un- 
der relative rotation between the real and spin space. 
Accordingly in below the spin quantization axis will be 
chosen for the most convenient presentation (and always 
along the total magnetization if it is finite). In partic- 
ular the configurations presented below do not rely on 
any special relative orientation between the net magne- 
tization and the rotational axis (which is always chosen 



a,s z). I shall also discuss the local magnetization den- 
sity m. The projection of m along a general direction 
is measureable in BEC experiments since it is given by 
the difference in the number density between the u and d 
species when one uses that direction as the quantization 
axis. 

Setting the variation of the free energy with respect 
to \&* to zero, one obtains the familiar Gross-Pitaevskii 
(GP) equations (generalized due to the presence of mul- 
tiple spin species P,f|). If c 2 > such as in the case 
of 23 Na, in the absence of a net magnetization the or- 
der parameter can be chosen so that only \^o is finite 
and obeys the GP equation in the usual form: = 
(-2177 V 2 + V ~ /i )*° + c ol*o| 2 *o- The order param- 
eter profile for \&o would then be completely analogous 
to that of a scalar order parameter with the interaction 
parameter given by Co (note then v&o is independent of Ci 
). In particular in the absence of any circulation and if 
one ignores the gradient term, ( the Thomas Fermi (TF) 
approximation ||) |^o| 2 = ^—^i^ — V). The structure 
of a singly quantized vortex would also be exactly analo- 
gous to that of a scalar order parameter investigated by, 
e.g., Dodd et al ||. For the discussions below it is also 
convenient to re-consider the same CV with quantization 
axis rotated by tt/2 about a horizontal axis. In this basis 
the above CV becomes two coinciding vortices of the u 
and d components with \^f u \ — Y&d\ an d their nodes at 
the trap center. We shall see below that in general the 
CV is very different from the ones just discussed. 

Typically in the experiments the cloud is trapped by 
a potential harmonic in all three spatial directions. For 
numerical simplicity I shall instead consider a cloud sub- 
ject to a harmonic potential only in the x — y plane but of 
uniform density within thickness d along the rotational 
z axis. It is reasonable to assume that the results be- 
low will be qualitatively applicable to a pancake shaped 
cloud trapped by a three dimensional, axially symmetric 
harmonic potential if the radii of the clouds perpendic- 
ular to the rotational axis are comparable. Rather than 
varying /i and SI to obtain a fixed total number of par- 
ticles and angular momentum, I shall simply present the 
types of CV for fixed /x's and Si's. However, I shall con- 
tinue to use the total magnetization (rather than H) as 
an independent variable |2| . I shall eliminate /i in favor 
of R = (2^/MaW 2 ) 1 / 2 , where lo is the (angular) trap fre- 
quency. In the absence of vortices and under the TF ap- 
proximation, the radius of the cloud and the total number 
of particles are independent of the value of c 2 and given 
by R and N = yf^(jp-) 4 respectively, where A is the 
size of the harmonic oscillator ground state wavefunction 
(A = (H/MaUJo) 1 / 2 ). Here a = (a ° + 3 2a2) is an effective 
scattering length for the interaction parameter cq. All 
the CV presented below has N « N . (|] 

I shall introduce the parameter e = (Xq/R) 2 which 
measures the deviation from the TF (e — > 0) limit, e de- 



pends only weakly on N for given trap parameters. I shall 
express vE^ in units of y/fj,/c (correspondingly the num- 
ber density n and the magnetization density rh in [i/c ) 
distances in units of R, and total particle number and 
magnetization mtot m units of N . With this, all physi- 
cal results depend only on the dimensionless parameters 
e, SI = 0,/ui o , c-2 = cijc and m tot . Anticipating future 
experiments on other atoms I will not fix c 2 to that of 
23 Na (though confining myself to c 2 > ). As a concrete 
example I shall consider mainly e = 0.1, SI = 0.45, ITc ] 
with the corresponding phase diagram shown in Fig. [I . 
I shall comment on other values of the parameters as I 
proceed. 

We begin by considering mtot = 0. I shall present 
the CV in two ways, each related to some of the CV 
structures discussed for mtot ^ below. The structure of 
the CV with quantization axis chosen so that it resembles 
most closely a vortex of ^o alone is as shown in Fig. [|. 
However, instead of an empty core, it is energetically 
favorable for some of the particles to convert to u and 
d species and appear near the center of the trap. For 
the present parameters, \^ u \ — \$d\ and each has two 
nodes with unit circulation. In Fig@ the order parameter 
along the line (chosen as the x axis) going through these 
singularities was shown. Note that the CV has broken 
cylindrical symmetry. pd| 

Another useful way of presenting the above CV is to 
use quantization axis rotated by 7r/2 about a horizontal 
axis with respect to those above. In this basis only ty u 
and ^d are finite. Each has one node, displaced by equal 
but opposite distance from the trap center (Fig pi). It 
follows that the local magnetization density is finite and 
points along the (present) z-axis, being negative for x < 
and positive for x > 0. Notice that at the singularity for 
say the (i-component, since |\I/d| = and \^ u \ ^ 0, locally 
the condensate is actually in the axial but not polar state 
(even though mtot and H are zero) 

With the use of this quantization axis we can also un- 
derstand easily the reason for the present CV structure. 
Due to the presence of the trap potential a vortex has 
maximum kinetic energy if its node is located at the cen- 
ter of the trap. It is thus energetically favorable for the 
nodes of the u and d components to move away from the 
trap center. MM For the system to be at an energy mini- 
mum, they move opposite to each other, creating regions 
where \^f u \ ^ \^d\, eventually balanced by the desire of 
the condensate to remain in the polar state. In this pic- 
ture it is obvious that L/N < 1. As SI increases, the 
^ u ,d singularities move closer to the center of the trap 
and L/N increases, [e.g., for c 2 = 0.2, L/N = 0.87(0.92) 
at SI = 0.45(0.5)]. The separation between these sin- 
gularities, and the region where the local magnetization 
is non-zero, increase with decreasing c 2 [correspondingly 
L/N decreases: e.g., at SI = 0.45, L/N = 0.84(0.80) for 
5 2 = 0.1(0.05)] 

Now we are ready to consider the structure of the CV 



with finite total magnetization. I shall describe each re- 
gion of the phase diagram Fig. \u in turn. 

I: In this region the favorable configuration is similar 
to that of Fig || except for an increase (decrease) in the 
amplitude of \& u and ^d (not shown) . The local and to- 
tal magnetization of this CV are always collinear. One 
can understand this configuration by considering the en- 
ergy under the magnetic field H. The CV has an order 
parameter and hence a magnetic susceptibility which is 
anisotropic. For a given magnitude of the magnetization, 
the energy is minimum if the direction of mtot is along 
that of largest susceptibility. The quantization axes used 
in Fig |2j and y above correspond to the principal direc- 
tions of the susceptibility tensor. It is intuitively rea- 
sonable that the CV has larger susceptibility along the 
quantization axis of Fig || (c.f. J|). 

II: For larger m to t the vortex of the minority species 
d disappears. The CV is replaced by a vortex of the u 
species with a d core. ( Fig. [|.) This can be understood 
by recognizing that the effective chemical potential for 
the d species is given by jx — H. Increasing m to t requires 
increasing H, hence decreasing \x — H. Eventually the 
effective chemical potential is too low to overcome the 
necessary kinetic energy required for forming a circulat- 
ing d component. This picture is supported by the fact 
that the critical mtot needed for the I > II transition in- 
creases with $7. The — J7L term in the free energy favors 
an order parameter with finite circulation, thus at higher 
angular velocity a larger H and hence mtot is required for 
the transition, [e.g, at £2 = 0.2, the critical m-tot ~ 0.2 
for Cl = 0.45 here (Fig Q) whereas m to t ~ 0.4 for Cl = 0.5] 

III: This occurs at still larger mtot and only for small 
ci- In this region the CV has a u vortex with a core 
filled by the species (Fig o) . The minimum magnitude 
of m to t needed for this new C V increases with 62 . These 
features can be understood by considering again the ef- 
fective chemical potential for the and d spins which are 
/1 and /1 — H respectively. The spin species is more 
favored by H, but suffers a stronger repulsion (than the 
d species) from the majority u species due to the spin 
dependent interaction C2(> 0). Only at sufficiently small 
C2 and large m to t does this CV become favorable. 

In Fig. O we display the local magnetization density 
m of this CV at points on the x axis, defined so that 
the phase difference between the u and components 
vanishes for x > 0. Near the trap center m points mainly 
along the horizontal, turning towards z, the direction of 
net magnetization, only further away. The magnitude as 
well as the z-component of m depend only on the radial 
distance from the center of the trap. The azimuthal angle 
of m is the negative of that of the corresponding physical 
point in space. It is interesting to note that the presence 
of the CV may not be apparent if one examines only the 
particle number density n, in strong contrast to the case 
of a scalar condensate. M . 

IV: This is the most intriguing region. At very small 



£2 (and not too small e's) the CV has spontaneous (spin) 
symmetry breaking in the sense that it has a net magne- 
tization even when H = 0. The configuration is similar 
to that of Fig. g except now the numbers of spin-up 
and spin-down particles are no longer equal (see Fig]?]). 
This configuration is stable (i.e. the topology of the CV 
remains the same except for a re-adjustment of the am- 
plitudes of ^'s) so long as the total magnetization is close 
to that of the 'spontaneous' one. The corresponding local 
magnetization density is as shown in Fig @. The direction 
of m thus rotates from —x through z to x as one moves 
along the physical z-axis. The existence of this sponta- 
neous mtot means more of the local order parameter is 
axial like, and thus this state is possible only for suffi- 
ciently small C2- For larger m to t, this configuration gives 
way to that of a u- vortex and a d core (Fig Q) , which has 
a larger susceptibility, discussed earlier. 

In conclusion I have shown that the internal vortex 
structure of a spin-1 Bose condensate in a harmonic trap 
is much richer than that of a condensate with a scalar 
order parameter. I thank T.-L. Ho for his comments on 
the manuscript. 
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FIG. 1. The phase diagram for e = 0.1, Q = 0.45. 
I (diamonds): vortices in both u and d, configuration 
as in Fig H except for a possible adjustment of relative 
magnitudes of |^ u ,d|; II (triangles right): vortex in u with 
d core as in Fig m III (triangles down): vortex in u with 
mainly core as in Fig H; IV (triangles up): configuration 
as in Fig. M. Dotted lines are guides to the eye. 
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FIG. 2. CV for c 2 = 0.2, fi = 0.45. mtot = 0. |*o| 
full line; |*„| = |* d |, dashed. 
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FIG. 3. The same CV as Fig |, but with different 

quantization axis. \9 U \, dashed-dotted; \^d\, dashed. 
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FIG. 4. CV for c 2 = 0.2, fi = 0.45, and mtot = 0.414. 



\^u\, dashed-dotted; |\Pd|, dashed. 
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FIG. 5. CV for c 2 = 0.02, Q = 0.45, and m tot = 0.77. 
|*u|, dashed-dotted; |* |, full-line. 
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FIG. 6. Local magnetization for the CV of Fig 
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FIG. 7. CV for c 2 = 0.02, fi = 0.45, and m tot = 0.03 
(H = 0). [*„], dashed-dotted; j* |, full line; |* d |, 
dashed. 
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FIG. 8. Local magnetization and density for the CV 
of Fig |. 
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Figure 1: CV for Fig. 2 of text, c 2 = 0.2, = 0.45, m tot = 0; |* t 
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Figure 2: Magnetization of the above CV, only x component is finite 



'PsLu' 

0.75 

0.6 

0.45 

0.3 

fe 15 

1 

0.5 



-0.5 

-1 

-1.5 




Figure 3: Contour plot of CV for Fig 3 of text, c 2 = 0.2, = 0.45, m tot = 0; 
i.e. same vortex as Fig 1 here but with different quantization axis. Note now 
* z = 
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Figure 4: 2D contour plot for CV for c 2 = 0.02, = 0.45, and m tc 
(Fig 7 of text) 
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Figure 5: Magnetization for the above CV. 



